RECENT QCD RESULTS IN u - N 
DEEP-INELASTIC-SCATTERING AT CCFR/NUTEV 



JAEHOON YU 
(for the CCFR/NuTeV Collaborations) 
MS309, FNAL, P.O.Box 500, Batavia, 
IL 60510, USA 
E-mail: yu@fnal.gov 

We present recent QCD results in v — N scattering at the Fermilab CCFR/NuTeV 
experiments. We present the latest Next-to-Next-Leading order strong coupling 
constant, a s , extracted from Gross-Llewellyn-Smith sum rule. The value of a s 
from this measurement, at the mass of Z boson, is a% JVXO (M§) = 0.114+g' u 5|. 
We also present a preliminary result of the CCFR F2 at large-x. This measurement 
of F2 is compared to other experiments and various models that provide nuclear 
effects. We discuss the previous strange sea quark measurement from the CCFR 
experiment and the prospects for improvements of the measurement in the NuTeV 
experiment. 



1 Introduction 

Neutrino-nucleon (f-N) deep inelastic scattering (DIS) experiments provide a 
good testing field for Quantum Chromo-Dynamics (QCD), the theory of strong 
interactions. The v-N DIS experiments probe the structure of nucleons and 
provide an opportunity to test QCD evolutions and to extract the strong cou- 
pling constant, a s . They are complementary measurements to charged lepton 
DIS experiments of nucleon structure functions. The advantage of z^-N DIS 
measurements over the charged lepton experiments is that u-N experiments 
can measure both F2(x,Q 2 ) and xFs(x,Q 2 ) due to pure V-A nature. The 
v — N differential cross sections are written, in terms of structure functions: 



dxdy ir 



Mxy y 2 l + AM 2 x 2 /Q 2 \ {v) 
V 2E„ 2 l + R(x,Q 2 ) J 2 



(i) 

where R(x,Q 2 ) = glI^t, the ratio of longitudinal to transverse absorption 
cross sections. The structure functions F%(x, Q 2 ) and xF^(x, Q 2 ) are extracted 
from fitting Eq. ^to the measured differential cross sections. 

In this paper, we present a next-to-next leading order (NNLO) determi- 
nation of a s from the Gross-Llewellyn-Smith (GLS) sum rule and a large- x 
structure function measurement from the CCFR experiment. We also discuss 
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the measurements of strange sea quark distributions from the CCFR experi- 
ment and the expected improvements of this measurement in NuTeV. 

2 The Experiments 

CCFR (Columbia-Chicago-Fermilab-Rochester)/NuTeV experiment is a v — N 
DIS experiment at the Tevatron in Fermilab. The CCFR experiment used a 
broad momentum beam of mixed neutrinos (y^) and antineutrinos (y^) from 
decays of the secondary pions and kaons, resulting from interactions of 800GeV 
primary protons on a Beryllium-Oxide (BeO) target. The NuTeV experiment, 
the successor of CCFR, used separate beams of or F M during a given running 
period, using a Sign-Selected-Quadrupole- Train (SSQT) 111. 

The CCFR/NuTeV detector 13 consists of two major components : target 
calorimeter and muon spectrometer. The target calorimeter is an iron-liquid- 
scintillator sampling calorimeter, interspersed with drift chambers to provide 
track information of the muons resulting from charged-current (CC) interac- 
tions where a charged weak boson (W + or W~) is exchanged between the 
(T7 M ) and the parton. The calorimeter provides dense material in the path of 
v (77) to increase the rate of neutrino interactions. The hadron energy resolu- 
tion of the calorimeter is measured from the test beam and is found to be : 
<j/E Had = (0.847 ± 0.015) /y/E Had (GeV) + (0.30 ± 0.12)/ E Had . 

The muon spectrometer is located immediately downstream of the target 
calorimeter and consists of three toroidal magnets and five drift chamber sta- 
tions to provide accurate measurements of muon momenta. The momentum 
resolution of the spectrometer (<r/p M ) is approximately 10.1% and the angular 
resolution is 9^ = 0.3 + 60/p At (mrad). 

3 a s from Gross-Llewellyn-Smith Sum Rule 

Once the structure function xF^ is extracted, one can use the Gross-Llewellyn- 
Smith (GLS) sum ruletj, which states that J xFs(dx/x) is the total number 
of valence quarks in a nuclcon, up to QCD corrections, to extract the strong 
coupling constant, a s . Since in leading order (LO), the structure function xF% 
is xq — xq, the valence quark distributions, integrating xF^ over x yields total 
number of valence quarks, 3. 

Since the GLS sum rule is a fundamental prediction of QCD and the 
integral only depends on valence quark distributions, a s can be determined 
without being affected by less well known gluon distributions. Moreover, since 
there are sufficient number of measurements of xF% in a wide range of Q 2 , one 
can measure a s as a function of Q 2 in low Q 2 where the values of a s depends 
most on Aqcd- 
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Figure 1: XF3 vs x for four lowest Q 2 bins. The solid circles represent the CCFR XF3 data 
and the inverse triangle represent the data from other experimental measurements. 

With next-to-next-to- leading order QCD corrections E, the GLS integral 
takes the form : 

f xF 3 (x, Q 2 )— = 3 (l - ^ - a(n,)(^) a - b(n f )A 3 ) AHT (2) 

Jq X \ 7T 7T 7T / 

where the term AHT is the corrections from higher-twist effects. The highex.- 
twist correction term, AHT, is predicted to be significant in some models 13, 
while others □ predict negligibly small corrections. We take AHT as one half 
the largest model prediction with the associated error covering the full range 
(AHT — (0.15±0.15)/Q 2 ). 

In order to perform the integration of xF 3 in the entire ranges of a;, we use 
data from other iy-IiDIS experiments (WA59, WA25, SKAT, FNAL-E180, and 
BEBC-GargamelleQ), to cover large- a; ranges that are not covered by CCFR 
due to geometric and kinematic acceptances. The CCFR data has a minimum 
x of roughly x = 0.002Q 2 . To extrapolate below this kinematic limit, we fit 
XF3 to a power law (Ax B ) using all the data points in x < 0.1. 

Since the data points do not populate the region x > 0.5 finely, it is 
necessary to interpolate between the data points. Thus, we use the shapes of 
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the charged lepton DIS F2 in x > 0.5, because the shapes of xF% and F2 should 
be the same in this region of x due to negligible contribution, from sea quarks. 
The charged lepton F2 data are corrected for nuclear effectsla before obtaining 
the shapes. Systematic uncertainties due to extrapolation and interpolation 
take into account the possible variations of the models in the low- a; region and 
the resonance peaks in charged- lepton F% shapes in large- a; regions. 

Figure |l| shows xF% in the four lowest Q 2 bins as a function of x. The solid 
circles represent the experimental xF^ data, The solid lines represent the power 
law fit used for extrapolation of xF% outside the kinematic limit of CCFR. The 
shaded area in each plot shows the region of x where the extrapolations are 
used for the integration. 

The xF% integrals are estimated in six Q 2 bins. The results in each Q 2 
bin are fit to the NNLO pQCD function and higher-twist effect in Eq. ||. This 
procedure yields a best fit of A^L NLO = 165MeV. Evolving this result to the 
mass of the Z boson, Mz, in NNLO, this corresponds to the value of a s : 

a NNLO [M 2 z) = .H4+o : 005 (stot) +0.007 (syst) ± . 05(iWy) (3) 

where systematic uncertainty includes the error in the fit and theory error 
represents the uncertainty due to higher twist effects and higher order QCD 
contributions. This result is consistent with the value extracted from the 
CCFR structure function measurement!]. 

Evolving to the mean value of Q 2 {— 3GeV 2 ) for this analysis results 
in af NLO {3GeV 2 ) = 0.278 ± 0.035 ± 0.05+°;of. If the higher-twist effect is 
neglected, the value of af NLO (M 2 z ) becomes 0.118. 



4 Large- 2; Structure Functions 

In a simple parton model, the contributions of a single stationary nucleons 
vanish as x, the fractional parton momentum, approaches to 1. However, the 
interplays between nucleons in a nuclear environment would add a momentum 
to initial nucleon that interacts with the neutrino, causing x to be greater than 
1. This nuclear effect has been observed and measured in EMC experiment 
and is called the EMC effect El. 

There have been several models to describe this nuclear effect : 11 Fermi- 
gas model which adds Fermi momentumE3 and a quasi-deutron effectEil, 2) few 
nucleon correlation modelt2l and 3) the multi-quark cluster modeM, which add 
exponentially falling F2 distributions. 

In order to explore these various models, the CCFR experiment measured 
large-x cross sections in the kinematic regions x > 0.75 and Q 2 > 50GeV 2 . 
This measured cross section is compared to a LO Monte Cairo predictions 
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Figure 2: Number of events vs x for all events in Q 2 > 50GeV 2 . The histogram represent 
the data and the points represents the fit results. 

that includes various detector smearing effects and incorporates various models 
separately. Figure ^| shows the number of events as a function of the scaling 
variable x in the region_r > 0.75. The structure function F 2 (x, Q 2 ) is then fit 
to the functional formEd F 2 (x, Q 2 ) = (1 - x) a (b + cx + dx 2 + ex 3 )(Q 2 ) f+9X g 
for x < 0.75 and F 2 {x,Q 2 ) = F 2 (x = 0.75, Q 2 ^-^- - 75 ) for x > 0.75. The 
resulting value of the parameter s from this fit is s — 8.3±0.7(stat)±0.3(syst). 

We also investigated a possible Q 2 dependence of this parameter s in two 
Q 2 bins, 50 < Q 2 < 100 GeV 2 and 100 < Q 2 < 400 GeV 2 . The resulting 
values of the parameter s in these two kinematic regions are s — 7.4 ± 0.9 
for lower Q 2 bin and s — 8.7 ± 0.9 in the higher Q 2 bin, showing only a 
weak Q 2 dependences. Table compares this result with other experimental 
measurements and models, demonstrating that the FNC prediction and the 
SLAC El 33 measurement of the parameter s agree with our result, while the 

Table 1: Comparisons of the exponent s of CCFR with other predictions and experimental 
measurements. 



Experiment /Predictions 


s 


CCFR 


8.3 ±0.8 


FNC Prediction 


8-9 


Baldin et al. Prediction 


- 6 


SLAC E133 


7-8 


BCDMS 


16.5 ±0.5 



a a:2.5693, b=0.2739, c=3.0437, d=-5.5172, e=2.5790, f=-0.0303, and g=-0.2185 



BCDMS result differs significantly. 

From this analysis, we find that the models without nuclear effect do not 
describe our data well. Incorporating only Fermi gas and quasi-deutron model 
is not enough to describe the excess of the data over x > 1. The FNC model 
describes the CCFR large- a; behavior well with the exponential fall off param- 
eter of s = 8.3 ± 0.9. This measurement indicates that thc-jauclcon structure 
function F2 at x — 1 does not vanish but has a finite value ta of ~ 2 x 10~ 3 . 



5 Strange Sea Measurements 

In v — N DIS, the flavor changing weak CC interactions provide a clean sig- 
nature for a scattering off an s quark (s — > c), resulting in a pair of oppo- 
sitely signed muons in the final state through the reaction v^iv^) + ./V — > 
+ c(c) + X where c(c) subsqeuantly undergoes a semileptonic decay 
c(c) — * s(s) + + V^v^). In contrast, detecting s quark with NC in- 

teractions in a charged lepton DIS is convoluted with s production and frag- 
mentation of strange mesons, requiring good particle identifications. The nu- 
cleon structure function F2 from v — N DIS together with that measured from 
charged lepton DIS could also measure s quark distributions. In addition, the 
difference of the v — N structure function xF% between v and V measures the 
nucleon strange quark contents as well. 

CCFR experiment performed both LO and NLO analyses 0. Since 
a neutrino from the semileptonic decay of the charm is involved in the fi- 
nal state, the CCFR measured the distributions of visible physical quantities; 

2 

x vis , E ms = E Had + p* + p + n 2 , and z vis = En ^ +p 2 which is the fractional 
momentum taken by the charm meson fragmentation. The measured distri- 
butions were compared to Monte Carlo predictions which incorporate detector 
acceptance and resolutions. 

The Monte Carlo included the singlet (xqsi(x, /s 2 )), non-singlet (xqNs( x , M 2 ) = 
xqv (x, M 2 )), and gluon parton distribution functions from CCFR structure 
function analysis assuming: 

xqv(x,fi 2 ) = xuy (x, fi 2 ) + xdy (x, (J, ) (4) 

xdv(x, /1 2 ) = Ad(l — x)xuv{x. fi 2 ) (5) 

xu(x, /1 2 ) = xus{x, /j 2 ) (6) 

xd(x, /1 2 ) = xds(x, /j 2 ) (7) 

The above assumes symmetric sea quark distributions. We then parameterized 
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the strange sea distribution to: 



xs(x, ll 2 ) = A s (l — x) 



a 



xu(x, ll 2 ) + xd(x, LL 2 ) 



(8) 
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where, A s is defined in terms of level and shape parameters, k and a. The 
level parameter n is defined as : 



The LO prediction included a simple parton model which has tree level 
calculations for scattering of a W off an s quark, resulting in a charm quark final 
state. On the other hand, the NLO prediction included the ACOT formalismll3 
whose calculation is performed in a variable flavor MS scheme. Specifically the 
NLO prediction included the Born and gluon-fusion diagrams. In addition, 
since the charm quark is massive, the threshold effect was taken into account 
in the prediction via leading order slow rescaling formalism where the scaling 
variable x is replaced with £ defined as : 



where m c is the mass of the charm quark. 

We then fit the Monte Carlo to data distributions of x V i S , E V i S , and z V i S 
for k, a, B c (the charmed meson semi-leptonic branching ratio), m c , and e 
(fragmentation parameter), using the values of CKM martix elements | V c d | 
and | V cs |, taken from PDG@. 

Table |^ summarizes the results of this analysis, comparing various frag- 
mentation functions at LO and NLO. The result shows that the strange sea 
level parameter, k, is in qualitative agreement between LO and NLO analyses 
while the shape parameter, a, is the same as other sea quarks in NLO but 
softer in LO. The charm quark mass differs between NLO and LO, due pre- 
sumably to the fact that the parameter m c is not a pole mass but rather a 
theoretical parameter that absorbs the lack of higher corrections. 

Aside from the large statistical uncertainty, there were two major sources 
of systematic uncertainties in performing this measurement in CCFR. First and 
foremost one is the identification of the secondary muons resulting from the 
decay of charmed mesons. Since the neutrino beam in CCFR was a mixture of 
and Vfj,, the experiment relied on transverse momentum of the muons rela- 
tive to the neutrino beam direction, , to distinguish the prompt muons from 



K 



Jo [xs(x,y 2 ) + xs(x,/j, 2 )} 
J [xu(x, /J, 2 ) + xd(x, fj,' 2 )] 



(9) 




(10) 
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Table 2: Summary of the CCFR strange sea analysis. The first set of the errors in each item 
is statistical and the second set is the systematic uncertainties. 





Fragmentation 


K 


a 


B c 


m c (GeV) 


NLO 


Collins- Spillcr 
e = 0.81 
±0.14 


0.477 

+0.046 
-0.044 
+0.023 
-0.024 


-0.02 

+0.60 
-0.54 
+0.28 
-0.26 


0.1091 

+0.0082 
-0.0074 
+0.0063 
-0.0051 


1.70 
±0.17 

+0.09 
-0.08 


NLO 


Peterson 
e = 0.20 
±0.04 


0.468 

+0.061 
-0.046 
+0.024 
-0.025 


-0.05 

+0.46 
-0.47 
+0.28 
-0.26 


0.1047 
±0.0076 

+0.0065 
-0.0052 


1.69 
±0.16 

+0.12 
-0.10 


LO 


Peterson 
e = 0.20 
±0.04 


0.373 

+0.048 
-0.041 

±0.018 


2.50 

+0.60 
-0.55 
+0.36 
-0.25 


0.1050 
±0.007 
±0.005 


1.31 

+0.20 
-0.22 
+0.12 
-0.11 



secondary muons. The second systematic uncertainty is the LO slow reseat- 
ing formalism to take into account the charm threshold effect in theoretical 
predictions. 

The experimental systematic uncertainty due to the identification of the 
secondary muon almost completely disappears in the NuTeV experiment due 
to the use of SSQT whose wrong sign neutrino contamination at a given mode 
is less than 10 -3 . In addition, the NuTeV experiment is working on improv- 
ing low energy muon energy measurements in the calorimeter using the in- 
tensive calibration beam program. In the theoretical predictions, NuTeV is 
planning to incorporate more complete NLO calculations and will investigate 
the dependence of the predictions on factorization scheme and various parton 
distribution function parameterizations. 

While statistical uncertainty dominated in the previous measurement, im- 
proving the systematic uncertainty would reduce the total error significantly. 
We expect the statistical uncertainty in this analysis would reduce by about 
a factor of 2.2 combining CCFR and NuTeV data. The NuTeV is currently 
planning to present a LO analysis by this summer. 

6 Conclusions 

CCFR measured GLS integral in six Q 2 bins and extracted the value of NNLO 
a s at the mass of the Z boson. The measured value of a s is: 

a™(Mf) = 1.114™ (11) 

which is consistent with the value measured from the structure function anal- 
ysis a and with world average. CCFR measured the large- x cross sections in 
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the region x > 0.75. Comparisons with various nuclear model seems to point 
that the FNC model along with Fermi-gas and quasi-deutron model describe 
the data well with the exponential fall off parameter s — 8.3 ± 0.9. The value 
of F 2 at x = 1 is measured to be F% Fe (x = 1) = 2.1 x 10" 3 . 

NuTeV has finished data taking in September 1997, running with separate 
beams of and , and a precision calibration (0.3%) is underway. Strange sea 
measurement at the NuTeV experiment is making progress, taking advantage 
of the separate beam and the intensive calibration. NuTeV expects a LO result 
in the strange sea measurement by the summer. 
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